Crossover of conductance and local density of states in a single-channel disordered 

quantum wire 
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The probability distribution of the mesoscopic local density of states (LDOS) for a single-channel 
disordered quantum wire with chiral symmetry is computed in two different geometries. An approxi- 
mate ansatz is proposed to describe the crossover of the probability distributions for the conductance 
and LDOS between the chiral and standard symmetry classes of a single- channel disordered quan- 
tum wire. The accuracy of this ansatz is discussed by comparison with a large-deviation ansatz 
introduced by H. Schomerus and M. Titov in Phys. Rev. B 67, 100201 (R) (2003). 
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I. INTRODUCTION 

The metal-insulator transition induced by disorder in 
the problem of Anderson localization remains poorly un- 
derstood from a theoretical point of view. On the one 
hand, it is true that the prediction for the existence of a 
metal-insulator transition induced by disorder has been 
verified by numerical simulations. On the other hand, 
the critical exponents at the metal-insulating transition 
measured numerically still cannot be extracted from the 
field theories which are believed to embody the metal- 
insulator transition^ The paradigm of this unpleasant 
situation is the plateau transition in the lowest Landau 
level of the integer quantum Hall effectual The develop- 
ment of reliable methods to compute analytically critical 
properties at a disorder-induced metal-insulator transi- 
tion remains to this date an open problem. 

The simplest known example of disorder-induced criti- 
cality in the context of Anderson localization was solved 
by Dyson in 1953i^ According to Dyson, the global den- 
sity of states (DOS) in the thermodynamic limit and close 
to the band center is anomalous when a single quan- 
tum particle hops with a random amplitude between the 
nearest- neighbor sites of a chain excluding any other form 
of disorder, in short, the one- dimensional random hop- 
ping problem. Subsequent works showed that Dyson's 
spectral anomaly is related to a diverging localization 
length )^!^!^ an anomalous decay of the envelope of wave 
functions^ and an anomalous probability distribution of 
the resistance upon approaching the band center^ 

The coalescence of the concept of universality imported 
from the theory of critical phenomena with the symmetry 
classification imported from random matrix theory has 
led to a classification of diffusive regimes and disorder- 
induced metal-insulator transitions in terms of 10 uni- 



versality classes that are uniquely characterized by the 
intrinsic symmetries preserved by the disorder (assumed 
weak) for any given dimensionality of spaced According 
to this classification, Dyson's singular DOS is a trade- 
mark of the chiral universality classes. The chiral uni- 
versality classes arc realized in problems of Anderson lo- 
calization for which the Hamiltonian anticommutes with 
some unitary operator for all realizations of the disorder 
and in which case the DOS is known to be singular at 
the band center in zero^iSiii quasi-one and two 
dimensions pi^iiStiiii^ Even though the chiral universality 
classes are critical at the band center, they have been far 
less studied than the orthogonal, unitary, and symplectic 
universality classes, which we will refer to as the stan- 
dard universality classes. For example, the probability 
distribution of the so-called mesoscopic local density of 
states (LDOS) has been known since 1989 from the work 
of Altshuler and Prigodin for a single chain with weak 
on-site disorder yiSiSLSi whereas it has not yet been com- 
puted for the one-dimensional random-hopping problem. 
The first aim of this paper is to fill this gap. 

If a diffusive regime exists, pcrturbativc techniques can 
be used to describe the crossovers between universality 
classesiSSiSi However, little is known quantitatively on 
the crossover between two different universality classes, 
here defined in terms of nonlinear-sigma models, say, 
from the diffusive to the localized regimes. This is even 
true in one dimension, which is by far the most stud- 
ied laboratory for Anderson localization as was recently 
illustrated by a flurry of works on the validity of one- 
parameter scaling in one dimensioni2^i2S*SSiSL2£ The ab- 
sence of a diffusive regime renders the concept of univer- 
sality classes ambiguous in one dimension in that sym- 
metry alone does not specify a universality class in one 
dimension. For example, the ratio of the localization 
length to the mean free path in a wire of finite width 



interpolates smoothly between its two limiting values for 
fully preserved or completely broken time-reversal sym- 
metry, respectivley, as a function of a weak magnetic 
fields Similarly, fine tuning of microscopic parameters 
is required to achieve delocalization (i.e., diverging lo- 
calization length) of quasiparticles at the fermi energy 
in a dirty superconducting wire of finite thickness with 
both broken spin-rotation and time-reversal symmetries 
whereas delocalization becomes generic in the thick wire 
scaling limit One can nevertheless define a symme- 
try class in one dimension by demanding that single- 
parameter scaling holds. Lack of universality in one di- 
mension can then be understood as the fact that there 
is no scaling limit for which single-parameter scaling be- 
comes a generic property of an ensemble of random mi- 
croscopic Hamiltonians of a given symmetry. Symme- 
try classes in one dimension cannot be construed as an 
enumeration of stable or unstable fixed points of some 
putative effective field theory. Computing the crossover 
between symmetry classes in one dimension nevertheless 
remains a well-defined problem. 

Time-reversal symmetry cannot be broken for a spin- 
less particle constrained to move in a onc-dimensional 
and simply connected world, in which case the only pos- 
sible crossover between symmetry classes takes place be- 
tween the chiral and standard classes. The second aim of 
this paper is to compute the crossover between these two 
symmetry classes in one dimension for the probability 
distributions of the conductance and of the LDOS from 
the ballistic to the localized regimes assuming a weak 
disorder at the microscopic level. In that regard, our re- 
sults are complementary to the analysis made by Schome- 
rus and Titov of the validity of one-parameter scaling in 
one dimension in which they computed exactly the first 
four cumulants of the logarithm of the conductance to 
leading order in the ratio of the length of the chain to 
the localization length assuming a weak disorder at the 
microscopic levelS*2& We thus propose an approximate 
ansatz for the crossover whereby the approximation con- 
sists of assuming that the phase and amplitude of the 
reflection coefficient separate. By comparing the cumu- 
lants of the Lyapunov exponent of the transfer matrix 
computed with our approximate ansatz to the ones com- 
puted with the large-deviation ansatz of Schomerus and 
Titov deep in the localization regime, we deduce that 
our ansatz captures the first cumulant, but fails with the 
higher cumulants. 

The paper is organized as follows. We formulate the 
one-dimcnsional problem of Anderson localization as a 
continuous Schrodinger equation with a relativistic ki- 
netic energy in Sec. ^1 A functional renormalization 
group equation is derived for the probability distribution 
of the refiection coefficient in Sec. IIIII Exact solutions 
to these coupled functional renormalization group equa- 
tions are given for the chiral and the standard classes, 
as well as an approximate solution that describes the 
crossovers between these two classes in Sec. IIVI The 
accuracy of the approximate solution for the crossover 
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FIG. 1: In this paper, we will consider two different simply- 
connected geometries for a strictly ID disordered quantum 
wire which are imposed by suitable boundary conditions at 
y — —L/2 and y — +L/2, respectively. In geometry (a) and 
(b) the disordered quantum wire is closed to the right and 
open to the left where it is connected to an ideal lead. The 
particle in geometry (b) is subjected to absorption whereas 
it is not in geometry (a). In geometry (c) the disordered 
quantum wire is open to the left and to the right where it is 
connected to ideal leads. 

regime is estimated in Sec. The probability distribu- 
tion of the LDOS is computed exactly in the chiral and 
standard classes whereas it is computed approximately 
in the crossover regime between these two classes in Sec. 
We conclude in Sec. |Vll| 



II. MODEL 

In this paper we want to investigate the statistical 
properties of the LDOS and dimensionless conductance 
g for a single quantum particle, which is (i) restricted to 
a simply connected and strictly one-dimensional (ID) ge- 
ometry and (ii) subjected to a static and weak random 
environment (disorder) that interpolates smoothly be- 
tween the standard and chiral symmetry classes. Condi- 
tion (i) implies that time-reversal symmetry is preserved 
for every realization of the disorder. In ID and assuming 
a metalfic ground state in the absence of the disorder, 
the ground state consists of two distinct Fermi points. 
Plane waves at the Fermi points are called left- and right- 
movers, respectively. In the standard symmetry class, the 
channels for forward and backward scatterings induced 
by an on-site disorder potential in the basis of left- and 
right-movers arc equally likely up to nonrandom oscilla- 
tory factors. In the chiral symmetry class, the disorder 
potential is off-diagonal in the basis of left- and right- 
movers. Presuming weak disorder in condition (ii) allows 
us to choose a kinetic energy that is a first-order differ- 
ential operator, for it only makes sense to linearize the 
spectrum about the two Fermi points when the disorder 
potential is weak. We thus choose to model the dynam- 
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ics of the single quantum particle and the static random 
environment by the Hamiltonian 



^ fj.=0 



(2.1a) 



where we have set the Fermi velocity v-p and H to one. 
The boundary conditions that will be imposed in the se- 
quel are depicted in Fig.^ They all obey condition (i). 
We have denoted by CTq the unit 2x2 matrix and by cT]^ 2 3 
the 2x2 Pauli matrices. The static random environment 
is represented by three independent potentials v^{y) e M, 
yLt = 0, 1, 2 that we choose, for mathematical convenience, 
to be supported on the interval —L/2 < y < +L/2, and 
to be white-noise and Gaussian distributed with vanish- 
ing mean, (disorder averaging is denoted by (...)) 

(y^iy)) = 0, (2.1b) 

{v^{yK{y'))^2g^S^J{y^y'). (2.1c) 

The three variances i 2 carry the dimensions of inverse 
length and need not be equal. A finite value 



My)) = A 



(2.2) 



can also be easily accommodated in our formalismt^ In 
this case, the clean system is insulating as A opens up 
a dimerization gap at the Fermi energy and the random 
environment mimics to a first approximation the static 
fluctuations of the phonons responsible for the dimer- 
ization sufficiently close to the Peierls transitioniS^ The 
problem of Anderson localization in the continuum de- 
fined by Eq. (|2.1() is a coarse-grained version of a single 
particle hopping on a simply-connected chain with a uni- 
form nearest-neighbor hopping amplitude subjected to 
weak and independent random fluctuations of an on-site 
potential and of the nearest-neighbor hopping amplitude. 

The generic symmetry obeyed by Eq. H2.1a|l is time- 
reversal invariance. 



H. 



If the condition 



.9o = .9i = 



(2.3) 



(2.4) 



is satisfied, or, for that matter, any equivalent condition 
obtained by a redefinition of the Pauli matrices through 
a position independent SU(2) rotation around (T3, Eq. 
(|2.1a|) has an additional chiral symmetry in that it anti- 
commutes with a^, 



a^Ha^ = -H. 



(2.5) 



Consequently, all nonvanishing eigenvalues of H come in 
pairs with opposite signs. The origin of the chiral sym- 
metry of the continuum model for any given realization 
of the disorder is a sublattice symmetry of an appropri- 
ate lattice regularization. A microscopic random Hamil- 
tonian realizing the chiral symmetry is, for example, a 



tight-binding Hamiltonian on a hypercubic lattice with 
strictly vanishing on-site energies and random hopping 
matrix elements restricted to nearest-neighbor sites, in 
short, the random hopping problem. 

The LDOS, DOS, conductance, shot-noise power, etc., 
can all be extracted from the elements of the 2x2 scat- 
tering matrix at the energy E that connects incoming 
("i") and outgoing ("o") states 



= S 



E 



i.K 



r t' 
t r' 



i.L 



(2.6) 



The amplitude for incoming and outgoing waves to the 
left- and right-hand side of the disordered region is de- 
noted by ■0'/°''"/^ here. The matrix elements r, r' G C 
are the reflection coefficients whereas the matrix elements 
t,t' £ C are the transmission coefficients. Alternatively, 
one can work with the 2x2 transfer matrix at the 
energy E which follows from the scattering matrix H2.6|) 
through the basis transformation implied by 



= Mi 



(2.7) 



Conservation of probability dictates that the scattering 
matrix is unitary or, equivalently, that the transfer ma- 
trix is pseudo-unitary 

{Se)^ a,SE = SEao{SE)^ =a„ (2.8a) 
(Me)^ ct,Me=Me<T3{Me)^ =a,. (2.8b) 

Conservation of probability also implies the existence of 
the polar decompositions 



Me = 



v'* 
u 



u 



' tanh X sech x 



„ , sech X tanh x I ^ \ u ' 

rj \ /is 



V 



(2.9a) 



cosh X sinh x 



u' I sinh X cosh x J ^ yO v' ^ ^ 



V 



(2.9b) 

where it^, u'e, Ve, v'e are independent complex numbers 



with 



1 and Xp £ 



The polar decomposition is not unique. For example, the 
sign of Xe can be absorbed into a redefinition of . 



. The variable 



'E' 

when restricted to the halfline 
[0, oo[, has the geometrical interpretation of a radial coor- 
dinate on a Riemannian manifoldi^^*^ The time-reversal 
symmetry of the Hamiltonian implies the transformation 
laws 



Se — (^e) ' 
a^M*E(T-^ = Me, 

which enforce the constraints u'p 



(2.10a) 
(2.10b) 



and v'e 



or, equivalently, = t^. The chiral symmetry of the 
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Hamiltonian implies the transformation laws 



-E) 



(2.11a) 
(2.11b) 



-E 



u 



_E and v'j^E 

'''*'E' = ''-*E' ^+E 



which enforce the constraints u'^ 
v_E or, equivalently, r_^E 

t'*E- Then, at the band center i? = with the condition 
(|2.4|) . the scattering matrix becomes Hermitian. 

Under the assumptions that (i) the disorder is weak 
and (ii) m^, u'e, Ve, v'e are all independently and uni- 
formly distributed on the unit circle in the complex plane, 
it was essentially shown in Ref. that the probability 
distribution A" of > obeys the Fokker-Planck equa- 
tion 

dX{x;L) 



dL 



Id d 

— — sinh(2a;)— csch(2a;)A'(a;;L) (2.12) 
M ox ox 



in the geometry of Fig. The mean free path £ is 
some function of ^g, gi, g2, and e. In a more gen- 
eral context of multi channel quantum wires assumption 
(ii) is known as the isotropy assumptions^ whereas the 
Fokker-Planck equation 12.12(1 is a special case of the 
Dorokhov-Mello-Pichard-Kumar (DMPK) equationi2Si2i 
A consequence of assumption (ii) is one-parameter scal- 
ing as encoded by the Fokker-Planck equation H2.12|l . 
i.e., the probability distribution of x and, consequently, 
of g = sech^ x depends on the single dimensionless pa- 
rameter L/i. It is evident that assumption (ii) breaks 
down at the band center E = and with the chiral con- 
dition 1)2.4(1 because Ue^q = w^^g and Ve^q = v'e^q 
are then both real valued. Correspondingly, the Fokker- 
Planck equation obeyed by the probability distribution 
X of — cxD < Xe=o < +00 is different from Eq. 1(2.12(1 and 
given by the diffusion equation 

dX{x;L) 1 d'^X{x;L) 



dL 



2£ dx^ 



(2.13) 



in the geometry of Fig.^iSi^i*2£ The mean free path £ is 
some function of 52- Again, the Fokker-Planck equation 
((2.13(1 encodes one-parameter scaling as the probability 
distribution of x and, consequently, of g = sech^ x de- 
pends on the single dimensionless parameter L/£. There 
arc no known multiple-parameter scaling equations for 
the LDOS, DOS, conductance, etc., which describe the 
crossover between the two limiting cases described by 
Eqs. 1(2.12(1 and 1(2.13(1 for values of L ranging from the 
ballistic {L smaller than the mean free path) to the lo- 
calized regime {L larger than the mean free path) to the 
best of our knowledge (see Refs. IsflEcllill for discussions 
of two-parameter scaling) . We fill this gap in the remain- 
der of the paper. 
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FIG. 2: (a) A thin slice of length 5L with a SL <^ £ ^ L 
is added to the left of the disordered region of length L. (b) 
Two disordered regions 1 and 2 with scattering matrices Si 
and 5*2, respectively, in a quantum wire. 



the LDOS v follow from the knowledge of the statistical 
properties of the reflection coefficient 



R exp {iq 



(3.1) 



which has been decomposed into its square modulus < 
R < 1 and phase < (j) < 2tt. For a disordered wire 
of length L + SL as depicted in Fig. |21 the reflection 
coefficient r^^sL related to the entries of the scattering 
matrix for the slice of length SL and to the entry 
of the scattering matrix Sj^ for the much longer segment 
of length L by the composition law 



' L+SL 



SL- 



(3.2) 



When the width 5L of the slice is much larger than the 
lattice spacing a but much smaller than the mean free 
path £, here defined bji2^ 



we infer the continuous Langevin process 



dR 



(3.3) 



— = -4£"i? + 2\/!R(l-i?)(t>isin0-W2Cos0),(3.4a) 

QjLj 



— = 2(f ' + wq) + (wi cosc/) + W2 sine 
dL 



by using the relations 



R 



1 



R, 
(3.4b) 



FUNCTIONAL RENORMALIZATION 
GROUP EQUATIONS 



III. 



We shall show in Sees. IIVI and IVII that the statisti- 
cal properties of the dimensionless conductance g and of 



= i^Vi - V2) SL, = l + i (I'o + E) SL, 

r'sE = (^Vi+V2)SL, t'sE = l + t{vo+E)SL, (3.5) 

which are valid up to first order in the disorder potentials 
Vg 2 J by expanding the right-hand side of Eq. 13.21) 
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to the same order. We note in passing that the mean free 
path is simply given by 

(3.6) 



2(51+52) 

in the Born approximation. For later eonvenienee, we 
have continued the energy E to the upper part of the 
complex plane, E = £' + i£" , £" > 0. The contin- 
uous Langevin process H3.4|l can also be formulated as 
the Fokker-Planck equation obeyed by the joint probabil- 
ity distribution function P(R, (j); L). Using the standard 
methods of chapter 3 in Ref. |43 say, one finds^ 



G 



dP 

^'t'dRd4) 



(3.7a) 



d^P 



J,2 ' 



where 



Fo(i?,, 
FJR,, 



and 



Grr{R, 



2uj + 2(2i? - 1) - 6Ci? cos20, 
-2ujR- {R-l){bR-l) 
-6Ci?(i?- I)cos20, (3.7b) 



2R{1 - Rf - 2Ci?(l - Rf cos 20, 
C(l - i?^)sin2(/), (3.7c) 



G..(i?,0) = 4Co + - U + 2 



2 

R + 2 



R 



1 

'r 

cos 2(f). 



We have introduced the short-hand notations 



9± 5i ±52: 
£ := ^7.9+, 
C := .9-/5+, 



t:^L/t, 

Co := %/g+- 



(3.8) 



Observe that the dependence on the polar angle arises 
solely from cos 20 and sin 20 in the Fokker-Planck equa- 
tion (|3.7|l . Thus, provided initial and boundary condi- 
tions are also periodic on the interval [0, 7r[, the joint 
probability distribution P(i?, 0; L) is periodic on the in- 
terval [0, 7r[ even though the phase is originally de- 
fined modulo 2tt. The periodicity tt is a consequence of 
the Langevin process (|3.4|) being invariant under Vi? 
— Vi? and ^ + TT. A dimerization gap (^2(2/)) = A 
produces the changes 

F^{R,(t>) ^ Fo{R,(t)) + — (i?,- I)^cos0, (3.9a) 
5+ V R 

Fj^{R,4)) Fj^{R,(t>) + — {R-\)^/R cos, 4>, (3.9b) 
9+ 

F^{R, cj,) F^{R, cjy) + A(i? + 1)1 sin0. (3.9c) 
^ 2g^ y/R 



The polar decomposition (|2.9|) suggests the change of 
variable 



^i? = tanh|a;|, (3.10) 
in terms of which the Langevin process H3.4|) becomes 
dx 



dL 



= —£ sinh 2x + vi sin - 



W2 COS ( 



(wi COS + W2 sin ( 



(3.11a) 



(3.11b) 



In turn, the joint probability distribution function 
W{x, 0; L) obeys the Fokker-Planck equation 



dW 
~dt 



J,W-J, 



dW dW 



(3.12a) 



where 
and 



1 

2 dx^ 



20-1 cosh 2x 
1 



x4> 



d^W 1 
d^^2 



1 



< 



sinh 2x 
uj sinh 2x + 



sinh 2x 



tanh 2x 
1 1 



cos 20, (3.12b) 



2 tanh 2x 2 tanh 2x 



cos 20, 



e + C 



a 



1 



c 



tanh 2x 
4Co ' 



sinh 2x tanh 2x 

- ( cos 20) , 
sin 20, 



sin 20, 



+ c- 



tanh-" 2x tanh^ 2 



(3.12c) 
■ cos 20. 



The Langevin process (|3.11() or the Fokker-Planck equa- 
tion (|3.12|) are invariant under a; — s- — x and — > -|- tt. 

Equations (|3.7|) or H3.12|l are, from a conceptual point 
of view, the main result of this paper as they determine 
under what conditions one-parameter scaling holds in a 
ID weakly disordered quantum wire. For example, Eq. 
(|3.12|) encodes the functional renormalization group flow 
of the joint probability distribution of the radial coor- 
dinate |a;| = arctanh-\/i? and phase of the reflection 
coefficient r. If it is the property that the functional 
renormalization group flow of the probability distribu- 
tion of X depends solely on the dimensionless ratio L/£ 
aside from x, which is understood as one-parameter scal- 
ing for then two examples of one-parameter scaling 
can be constructed from Eq. H3.12|l as follows. First, the 
Fokker-Planck (DMPK) equation follows from in- 

sertion of the ansatz 



Wix,^;L) 



1 

2^ 



X{x;L) 



(3.13) 
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into Eq. (|3.12() with w = and from integration over 
(j) G [0, 2ir[ of the resulting equation. Second, the Fokker- 
Planck (DMPK) equation H2.13|l follows from insertion of 

the ansatz 



0) + 6{(f>-TT)]X{x;L) (3.14) 



with gQ = g^=uj~e~0 into Eq. H3.12|l and from inte- 
gration over (j) of the resulting equation. One-parameter 
scaling is not the rule for generic initial and boundary 
conditions of Eq. H3.12|l and for generic values of e and 
5o 1 2 due to the absence of a diffusive regime in IDm^ 
This is most clearly seen by the fact that the localization 
length ^, as defined by the typical dependence of the con- 
ductance g = sech^ x on L ^ £, is finite. ^ = 2€, for the 
standard symmetry class, but diverges in the chiral limit 
gQ ^ gi = Lu = e = 0. Hence, the localization length 
must depend on a second microscopic parameter aside 
from £ for generic values of £ and 5o i 2- Equations H3.7() 
and H3.12|l describe the full crossover from the chiral to 
the standard class with no restrictions on the values of x 
and L. 

For any given absorption lu > 0, the parameter space 
(C, Co> e) with Co > and C, £ e M of Eqs. and ^01^ 
is three-dimensional and is depicted in Fig. |21 Remark- 
ably, the parameter ( always enters Eqs. H3.7|l and (|3.12|l 
as a prefactor to cos 2(j) or sin 2(j) and conversely cos 2(f> or 
sin 2(j) are always multiplied by As we shall see in Sec. 
II VI this property turns out to be crucial in our study of 
the crossover between the chiral and standard symmetry 
classes. 

Equation (j^.Tf) was derived in the standard class by 
Abrikosov, Melnikov, and Kumar, and by Rammal and 
Doucot among others starting from the continuous non- 
relativistic Schrodinger equation4^i^^i^S*^ The large x 
limit 



dW 
dt 



t^e^^ - 2Ccos2<?!)) W 



^22- 



1 3, 
- + -Ccos2<?i 



dW 
dx 



dW 

+ (-e-3Csin2^) — 
+ -(1-Ccos20)^ 



-hCsin(2(/)) 



dxd(j) 



-h(2Co + l + Ccos2</>) 



(3.15) 



of Eq. (|3.12() was derived by Schomerus and Titov when 
= and t/g = g^ starting from a Langevin process for 
the eigenfunctions of the discrete Schrodinger equation 
on a chain 




FIG. 3: Parameter space for the Fokker-Planck equations 
13.711 or 13.121 1. The point Chi has the coordinates £ = ("q = 0, 
C = — 1. It represents the chiral symmetry class. The point 
KWl has the coordinates e = 0, Co = ~C = 1- H represents 
the Kappus-Wegner anomaly. The two shaded planes repre- 
sent the standard symmetry class. The regions C, > and 
C < are equivalent in that they are related by a rotation by 
an angle tt around 0-3 of the Pauli matrices in the Hamiltonian 
12.111 . Especially, this rotation relates Chi and KWl to the 
points Ch2 and KW2, respectively. 



IV. SOLUTIONS TO THE FUNCTIONAL 
RENORMALIZATION GROUP EQUATIONS 

Exact and approximate solutions to the Fokker-Planck 
equation H3.7|) are derived in this section. We will con- 
sider three regimes of parameter in Fig. |21 The first 
regime corresponds to the chiral symmetry class. The 
second regime corresponds to the standard symmetry 
class. The third regime describes a specific path con- 
necting the two symmetry classes. 

The chiral symmetry class is given by the condition 



0, 



(4.1a) 



which defines the band center of the energy eigenvalue 
spectrum of Hamiltonian (|2.1a() when 



c 



-1, 



Co = 0, (.go = .gi = 0), (4.1b) 



i.e., when Hamiltonian H2.1a(l anticommutcs with Ui as in 
Eq. (|2.5|) . The standard symmetry class occurs whenever 



or 



C = (.gi = 52) 
1/Co = 0, 

i/kl = o 



(4.2a) 
(4.2b) 
(4.2c) 
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since any one of these three conditions decouples R from 
(j) in Eq. Htj.7|l and is consistent with a stationary proba- 
biUty distribution of 4> which is uniformly distributed. In 
this paper we shall only consider the crossover between 
the chiral and standard symmetry classes along the line 



c 



-1, 



Co = 0, 



< e < oo. 



(4.3) 



We must naturally choose the initial and boundary 
conditions to be imposed on the solutions to the Fokker- 
Planck equation (|3.7|) . Our choice is motivated by the 
computation of the probability distribution of the LDOS 
that we will carry out in Sec. IVll Two geometries de- 
picted in Fig. El will be considered in Sec. IVII The rel- 
evant geometries for the probability distribution of the 
reflection coefficient r that enter the computation of the 
LDOS are depicted in Fig.Q In Figs. Hi and^) the disor- 
dered quantum wire is closed on the right-hand side and 
connected to a reservoir on the left-hand side. In Fig.^ 
conservation of probability holds. In Fig.^ conservation 
of probability does not hold due to the presence of an ab- 
sorption encoded by a finite imaginary part £" ~ Log^ of 
E. In Fig. ^ the disordered quantum wire is open at 
both ends and fiux conservation (no absorption) is as- 
sumed. We choose the ideal initial condition r = 1 and 
r ~ when L ~ for the semiopen and fully open ge- 
ometries of Fig. ^ respectively. In the chiral symmetry 
class, we demand that the phase of r is either or tt with 
probability 1/2 initially. In the standard symmetry class, 
we demand that the phase of r is uniformly distributed 
in the interval [0, 27r[ initially. In view of the periodic- 
ity of the Fokker-Planck equation (|3.7|) we always impose 
initial conditions with a periodicity of tt. 

We shall now construct exact solutions to the Fokker- 
Planck equation H3.7|l in the chiral and standard sym- 
metry classes, respectively, for the three geometries of 
Fig. n We shall also solve the Fokker-Planck equation 
(|3.7|) in the crossover regime (|4.3f) for the geometry of 
Fig. nji. We could not construct exact solutions to the 
Fokker-Planck equation (|3.7|) in the crossover between 
the chiral and standard symmetry classes for the geome- 
tries of Figs, and We shall construct, however, 
an approximate solution to the Fokker-Planck equation 
(|3.7|) in a crossover regime that interpolates between the 
exact solutions for the geometries of Figs.^^ and^J- The 
accuracy of this approximate solution will be discussed 
in Sec. El 



A. Disordered quantum wire closed on the 
right-hand side without absorption 

In geometry^, an incoming plane wave from an ideal 
lead is perfectly refiected by the disordered quantum 
wire. Hence, the reflection coefficient r of the disordered 
quantum wire must be a pure phase, rr* = 1, for all 
t = L/£. We thus insert the separation-of- variable ansatz 



into the Fokker-Planck equation H3.7() that, after integra- 
tion over R, reduces to 



5$ 5$ 

— = -2Ccos(20)$-[£ + 3Csin(2^)] — 



•[2Co + l + Ccos(20)] 



92$ 



(4.4b) 



By assumption, we have switched off the absorption cu = 
0. Our choice of initial condition shall depend on the 
proximity to the parameter regimes (HU and (I12J but 
shall always be periodic on the interval [0, tt[. 



1. Chiral symmetry class 

In the chiral symmetry class (|4.1f) the Fokker-Planck 
equation l)4.4b|l reduces to 

— = 2cos(2(/.)$ + 3sin(2(/.)— + [l-cos(20)]^ 

(4.5a) 

and has the normalized, stationary, and periodic solution 



1 



$(0) ^ ^(0 _ 0) + -(5((^ - tt) 



(4.5b) 



whenever the phase is or tt with probability 1/2 ini- 
tially. The same conclusion could have been anticipated 
by inspection of the Langevin process 13.4|l or from the 
fact that the chiral symmetry class (|4.1|l imposes the con- 
dition that the scattering matrix is Hermitian, see Eq. 



Standard symmetry class 



In the standard symmetry class H4.2a|l . cos 20 and 
sin2(/) drop out from the Fokker-Planck equation (|4.4b|l : 



5$ 

'dt 



+ (2Co + 1) 



92$ 



(4.6a) 



Its normalized, stationary, and periodic solution is 

1 



(4.6b) 



P{R,(t);t) = 5{R-l)^{(l)-t) 



(4.4a) 



whenever the phase ^ of r is uniformly distributed ini- 
tially. More generally, the solution to Eq. H4.6a|l con- 
verges to the uniform probability distribution (|4.6bl) as 
soon as t 3> 1 for any given initial probability distribu- 
tion. This is also true for the standard symmetry class 
(|4.2b|) and H4.2c|) as cos(2(/)) and sin(20) are ineffective in 
these limits. 

When ^ 52 ^-i^d |e| < 1 the stationary probabihty 
distribution of r need not be uniform in [0, 27r[. For exam- 
ple, consider the case C = +1, e = 0, Co = 5(1 + C) = 1j 
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for which the Fokker-Planck equation Ij4.4b|) reduces to 



(9$ a$ 

— = -2cos(20)$-3sin(20) — 
+ [3 + cos(20)] — 



(4.7a) 



and has the stationary and normahzed sohitionSSi^ 

\/2^ 



r2(i/4)yr 



cos^ 



(4.7b) 



with a periodicity of tt. Here T{z) denotes the vahie of 
the gamma function for z € C. Equation (|4.7b|l is a 
special case of the stationary solution 



$((^) oc 



1 



(4.7c) 



V2Co + l + Ccos(2(/.) 

to Eq. (|4.4bp with e = as shown by TitoviiS In the 
chiral limit ^ — 1 and Cq — > 0, this solution becomes 
unnormalizablc and we have to use Eq. Ij4.5b|) as a sta- 
tionary solution instead. 



3. 



Crossover regime 



For simplicity, we consider the crossover regime (|4.3|l 
The chiral symmetry class (|4.1|) is reached when e = 0. 
The standard symmetry class (|4.2c|l is reached when 
jej ^ 1. In the crossover regime, the Fokker-Planck equa- 
tion Il4.4bll reduces to 



$. (4.8) 



A stationary solution to Eq. H4.8() is constructed from the 
solution 



$(0) 



sm ( 



4> 



^+(£/2)(cot 0'— cot <; 



(4.9) 



valid for e > and on the interval [0, tt[, by periodic ex- 
tension to the interval [0, 2tt[. The constant M is chosen 
so that the solution 1)4. 9|l is properly normalized. The 
solution interpolates between Eqs. (|4.5b|) and l)4.6b|l as 
e ^ and e ^ 00, respectively, as it should be. Equa- 
tion (|4.9fl is an important intermediary result as we will 
use it to construct an approximation to the probabil- 
ity distributions for the conductance and LDOS in the 
crossover regime that interpolate between the exact chi- 
ral and standard limiting probability distributions. 



B. Disordered quantum wire closed on the 
right-hand side with absorption 

In geometry^, an incoming plane wave from an ideal 
lead is not perfectly reflected by the disordered quantum 
wire due to a finite absorption ui = £" jg^^. Hence, the 
reflection coefficient r of the disordered quantum wire is 
not a pure phase anymore, r = Vi?cxp(it/)). 



1. Chiral symmetry class 

In the chiral symmetry class (|4.1ll the phase <j> and 
squared magnitude R of the reflection coefficient r must 
necessarily separate since < r < 1 is necessarily real 
valued, 



P(i?, (I);t) 



7^(i?;^). (4.10) 



After insertion of Eq. 1)4.10)1 into the Fokker-Planck equa- 
tion ()3.7)l one finds that TZ{R;t) obeys 



dn 

'dt 



2uj + 2{3R-2) 



+2R{1 - RY 
with the stationary solution 



2ljR + 3{1- R){l-3R) 



dn 

'dR 



7^(i^) = A/"- 



exp 



R{1 - R) 



(4.11a) 



(4.11b) 



The constant J\f is chosen so that the solution ()4.11b)l 
is properly normalized (this is always possible for w > 
0). The stationary solution l)4.11bp has a square-root 
singularity at i? = 0, an essential singularity at i? = 1, 
and is normalizable for R £ [0,1]. We shall see below 
that the probability for R to be in the vicinity of is 
enhanced in the chiral symmetry class 1)4.1)) compared to 
the standard symmetry class (14.2)1 . 



2. Standard symmetry class 

In the standard symmetry class 1)4.2)1 the phase and 
squared magnitude R of the reflection coefficient r can 
reasonably be taken to separate in view of the initial 
condition of a uniformly distributed phase 



P(i?,(/);i) = ^niR;t). 



(4.12) 



After insertion of the separation-of- variable ansatz ()4.12)) 
into the Fokker-Planck equation ()3.7)l one flnds that 
TZ{R;t) obeys 



dn 

'dt 



2uj + 2{2R- 1) 



n 



2ljR+{1 - R){1 

d^n 



5i?) 



an 

'dR 



-R{1 - RY 



OR^ 



with the stationary and normalized solution 

cxp I 



n{R) = 2Lue^ 



2u) 



{l-Rf 



(4.13a) 



(4.13b) 
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The stationary solution (|4.13b|l has an essential singu- 
larity when R = 1. The stationary solution (|4.13b|) was 
found by Pradhan and Kumar in Ref. |^ When the ini- 
tial condition on the probability distribution of 4> is not 
that of a uniform distribution, it is presumed that the 
deviations of the solution to the Fokker-Planck equation 
(13. 7|) from the scparation-of-variablc ansatz H4.12|l with 
the stationary solution (j4.13b|) are short transients. We 
do not know of an analytical verification of this assump- 
tion. 



3. 



Crossover regime 



There is no compelling reason to believe that the 
squared amplitude R and the phase (j) of the reflection 
coefficient r separate in the crossover regime. However, 
short of an explicit exact solution to the Fokker-Planck 
equation 1)3. 7|l in the crossover regime, we shall never- 
theless pursue a strategy relying on an approximation 
built on a separation-of-variable ansatz that allows us to 
interpolate between the chiral and standard symmetry 
classes. One possible alternative to the separation-of- 
variable ansatz that we tried is a perturbative expansion 
about the standard symmetry class (|4.2|l using a Fourier 
expansion of the joint probability distribution for R and 
(j). However, this approach has the drawback that it can- 
not be expected to interpolate all the way to the chiral 
symmetry class (|4.1ll in view of the nature of the singular- 
ities characterizing the limiting probability distributions. 
For this reason wc will not present the perturbative ap- 
proach in this paper. 

Our starting point is the separation-of-variable ansatz 



P{R,^;t)=niR-,tMc^;t). 



(4.14a) 



We further assume that the probability distribution func- 
tion for the phase (f> shows a quick relaxation to its 
stationary solution, i.e., is assumed to be t- 

independent. 



(4.14b) 



Insertion of Eq. (|4.14|) into the Fokker-Planck equation 
(|3.7(l followed by an integration over (j> yields 



on 

'dt 



2lu + 2{2R - 1) - 2 Ca{R - 1) 11 
2loR+{1 - R){1 - 5R) 

- 2Ca(l - i?)(l - 2R) 



(1 - Ca)i?(l - RY 



dn 
'Wr 

(4.15a) 



where the real constant a is the ratio of two Fourier ex- 
pansion coefficients, 

„ ■f"-'?-'^^''^'^', (4,15b) 



A stationary solution to Eq. H4.15a|l is 



n{R) =N' 



exp 



(1-Ca)(l--R) 



^-C"/(i-C") [l-R) 



2/(l-Ca)' 



(4.15c) 



This solution reduces to the known results in the stan- 
dard [Eq. Ij4.13b|l ] and chiral [Eq. (j4.11bl) ] symmetry 
classes. In the standard symmetry class ()4.2|l the prob- 
ability distribution of the phase is uniform <&(</)) = l/27r 
and hence C^a = 0, whereas in the chiral symmetry class 
(14.1(1 it is a sum of the two delta functions <&(0) ~ 
i(S(</) - Q) + \5{<i) - tt), a = 1, C = -1- There re- 
mains considerable arbitrariness in the choice of $ be- 
cause many different $ share the same a. In practice 
we will choose the stationary solution H4.9|) since wc then 
recover the true solution in the semiopcn geometry of 
Fig. ^ as w ^ 0. We postpone to Sec. the discus- 
sion of the accuracy of the approximation implied by the 
separation-of-variable ansatz ((4.15|l with the choice H4.9|l 
for the stationary probability distribution of 4>. 



C. Disordered quantum wire opened at both ends 
without absorption 

In geometry^, an incoming plane wave from an ideal 
lead is partially reflected and partially transmitted by the 
disordered quantum wire. Following Abrikosov in Ref. 
we trade the squared magnitude R of the reflection 
coefficient r for the resistance (inverse of conductance) 

Q := -r—^ e [l,cx)] (4.16a) 



l-R 



under which 



Q{Q,^-t)=P{R,cj)-t)^ 
ag 

obeys the Fokker-Planck equation 
dt 



(4.16b) 



(4.17a) 



where 



^(e,0) = 2c^(2e-l), 

A^{e,4>) = -2uJB{Q-l)-{'iQ-l), (4.17b) 

C2^?2-2£< + l . 



2 q{q-1) 



■ sin 2(j), 



and 



b,,{bA) = 2^?(^,-l)(l-Ccos20), 

B,4q,<I>) = C(2f?-l)sin2,^, 
1(2^-1)2 



(4.17c) 



2 gig-1) 
C(2g-1)2 



2 gig-1) 



cos 20. 
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Chiral symmetry class 



2. Standard symmetry class 



In the chiral symmetry class (|4.1|l insertion of the 
separation-of-variable ansatz 



In the standard symmetry class 1)4. 2|) insertion of the 
separation-of-variablc ansatz 



Q{Q,(t>;i) = 



i<5(^-0) + i<5(^-7r) 



Q{Q-t) (4.18) 



into the Fokker-Planck equation H4.17|l yields 

= 2Q + 3(2,-l)^ + 2,(,-l)^ 



dt 

with the initial condition 

Q{Q-,t^Q) = 5{Q-l). 



(4.19a) 



(4.19b) 



The normalized solution to Eq. H4.19|l has the integral 
representation 



1 «-i/2 



du 



TT \/47rt'3 J a (1 - u)^Q - u 



^/v{l - v) 



with the auxiliary function and variable 

Q — u 



(4.20a) 



D(u, V] g) :~ In 



vu{l — u) 



t' := 2t. (4.20b) 



The integrations over u and v can be performed in closed 
form, 



exp 



(arccosh g^^'^ 
2t 



(4.20c) 



By shifting the minimum of the resistance g from 1 to 0, 
p := g—l,we recover the same asymptotics for the proba- 
bility distribution of p computed by Stone and Joanopou- 
los in Ref.|2 with a log- normal distribution of the nearest- 
neighbor hopping amplitudes of the discrete Schrodinger 
equation when p — > or p — > oo. The change of variable 
g ~ 1/ g yields the probability distribution 



1 



exp 



(arccosh 5 



-l/2^ 



2t 



(4.21) 



for the conductance g in the chiral symmetry class 
(|4.1|) . The same result follows from solving the diffusion 
(DMPK) equation H2.13|l and performing the change of 
variables g = 1/ cosh^ xi^ 



into the Fokker-Planck equation H4.17(l yields 



(4.22) 



dQ dQ d'^Q 

^ = (2.-l)^ + .(.-l)^ (4.23a) 



with the initial condition 

Q{g-t = Q)^5{g-l). 



(4.23b) 



As shown by Abrikosov in Ref. 143. the normalized solu- 
tion to Eq. H4.23|l has the integral representation 



-t/A 



du ■ 



-[r>o(«;e)]V4t 



VIttF Jo ["(1 - u 



— --7^Do(u;£>) 




ye 



dy 

( cosh^ y - g) 



with y 



iln 



2g (l - V^-Q-') 
Dq{u] g) = In 



1/2 



and 



(4.24) 



g — u 



u{l — u) 



(4.25) 



The change of variable g = 1/g yields the probability 
distribution 



5i9;t) = 



T^t^g^ Jy 



dy ■ 



ye 



-(t/4)-feV*) 



1 2 -,\l/2 

g cosh y — I) 



(4.26) 



with := — i In [2^ ^ (l — VI ~ s) ~ l] ^^i the conduc- 
tance g in the standard symmetry class H4.2|l . 



3. Crossover regime 



As was already the case in Sec. IIV B 31 we could not 
solve the Fokker-Planck equation (|4.17|) exactly in the 
crossover regime. We thus try the separation-of-variable 
ansatz 



(4.27) 



We are again assuming that the probability distribution 
$ is t independent and we made the specific choice for 
$ given by Eq. (|4.9() . Insertion of Eq. (|4.27() into the 
Fokker-Planck equation (|4.17|) followed by an integration 
over the phase (f) yields 



dQ dQ 

-g = -2CaQ+i2g^m-2Ca)^ 



d^Q 
dg^ 



(4.28a) 
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with the initial condition 

Qig;t = 0) = Sig-l). 



(4.28b) 



The normalized solution to Eq. H4.28|l has the integral 
representation 



sm 



—TrC,a \ g' 



„Ca/(l-Ca) 



(1 - l)Ca/{l-Ca} 

^ u I 



exp 



4(1 -Ca)2 



X I du 

lo [u{l-u){g 

»i 

X I dv 
lo 



u)]l/2 
-D^{u,v;g)/4t' 



1,1/2(1 _ yy/a-ca) 



D{u,v]g) 



(4.29a) 



with the auxiliary function and variable 



D{u,v] g) 



In 



vu{l — u) 



t' := (1 - Ca)t. (4.29b) 



As it should be, we recover Eq. (|4.20() when = —I 
and Eq. H4.24|l when a = 0. The approximate probabil- 
ity distribution Q for the conductance g in the crossover 
regime is obtained from Eq. (|4.29|) through 



Gi9;t) = g-^Q{g;t), 



Qia) 



(4.30) 



The dependence on t for the mean and variance of the 
conductance in the crossover regime computed from Eq. 
(|4.3U|I are depicted in Fig.0] In the chiral symmetry class 
C = — l,a = -|-l, the mean and variance decay like t"-'^/^. 
For a = 0.5 and a = 0, the mean and variance decay 
exponentially with t. 



(a) 




(b) 



CO 

> 



0.10 



0.05 



0.00 




FIG. 4: The mean (a) and variance (b) of the conductance 
computed from the probability distribution (14.31)11 as a func- 
tion of t = L/t for a = 0, 0.5, 1 and C = -1- 



V. VALIDITY OF THE SEPARATION ANSATZ 
FOR THE CROSSOVER REGIME 

We have constructed an approximate solution to the 
Fokker-Planck equations H3.7|l and H4.17|l in the crossover 
regime (|4.3ll by assuming the separation of the squared 
magnitude R from the phase (j) of the reflection coefficient 
r. In this section we shall investigate the accuracy of this 
approximation to the crossover regime (|4.3(l for very large 
values of the length L of the disordered quantum wire, 
i.e., L much larger than the localization length ^, which 
will be defined below. To this end, we shall compare the 
cumulants of the radial coordinate x computed from the 
scparation-of-variable ansatz with the same cumulants 
computed from a large-deviation ansatz introduced by 
Schomcrus and Titov in Refs. l25l and l26l 



when t — !■ oo of the solutions to Eq. (|3.15|) in the crossover 
regime. We try the separation-of-variable ansatz 



W{x,d^;t)=X{x;t)<^{<j>), 



(5.1) 



where $ is the stationary and normalized solution (|4.9|1 . 
Insertion of Eq. H5.1|) into Eq. (|3.15|) with subsequent 
integration over <j) yields 



dX 
'dt 



1 / 

-2(1 + ^") 



1 



(1 - c«) 



dx"^ 



(5.2a) 



Equation H5.2a|l is nothing but the diffusion equation 
with a constant drift proportional to the amount of chiral 
symmetry breaking (1 + C,a), i.e., it is solved by 



Large t limit in the crossover regime with the 
separation-of-variable ansatz 



X{x;t) = 



exp 



[.T-(l+Ca)t/2]" 



Vtt (1 - Ca) t 



(5.2b) 



We begin from Eq. (|3.15|) without absorption and in 
the crossover regime. We seek the asymptotic behavior 



normalized to the real line a; € M. In the approximation 
(I5.2b|) , the only nonvanishing cumulants are the first and 
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second cumulants, 



{X) 

1 + Ca 



1 + Ca, 



1 - Ca 



(5.3a) 



(5.3b) 



Provided (^a ^ —1 (i.e., away from the chiral symmetry 
class), the random variable x becomes self- averaging in 
the thermodynamic limit. It then makes sense to identify 
the localization length ^ through 



oo. 



which, together with Eq. (|5.3a() . gives 



2£ 



(5.4) 



(5.5) 



In the standard symmetry class (|4.2|l a = and the local- 
ization length is twice the mean free path ^ = 2 x ^ as is 
well knowniS In the chiral symmetry class (|4.1|) (a = —1 
and the localization length defined by Eq. H5.4|l is not 
a finite number anymore but a random variable as the 
right-hand side of Eq. H5.3a|l vanishes whereas the right- 
hand side of Eq. (|5.3b() remains finite. For the crossover 
regime (|4.3|l . the constant a implicitly depends on e. The 
asymptotic behavior of a for small e > can be evaluated 
as, with the probability distribution (|4.9|l . 



-f-oo 



"■^ (a;2 + l)i/2 J- 



dy 



{y2 + l)l/2^ 



a 



2/. 
1 



+°° rl-r 1 
oo "■^ (x2 + l)3/2 



oo (1/2 + 1)1/2 



■ Ine, 



=e(y-x)/2 



(5.6) 



since a logarithmic divergence occurs both in the x and 
y integrals in the numerator as we take £ — > 0, while it 
occurs only in the y integral in the denominator around 
y = — oo. It is then natural to identify 



2£ 



(at C = -1 and for e < 1) (5.7) 



with the typical localization length^ (~ | In jej | x ^) when 
approaching the chiral symmetry class H4.1|l through the 
crossover regime (|4.3|l . We have also computed the de- 
pendence of a on e by evaluating Eq. H5.6|) numerically. 



In Fig. El 1/(1 — a) is plotted as a function of e in the 
crossover regime We observe that 1/(1 — a) asymp- 

totically behaves as —(1/2) In \e\ in the vicinity (e < 0.1) 
of the chiral symmetry class (|4.1ll . 

12 




FIG. 5: Numerical evaluation of 1/(1 — a) as a function of 
e in the crossover regime 14.31 . The integral a defined in Eq. 
M.lSbl l is evaluated numerically for the stationary distribu- 
tion of the phase H4.9|l . The line representing —(1/2) log |e| is 
a guide to the eyes. 



B. Large t limit in the crossover regime with the 
large-deviation ansatz 

Schomcrus and Titov in Refs. l25l andl2^ devised a sys- 
tematic method to compute all the cumulants of the con- 
ductance and the LDOS to leading order in ^/L(^ 1). 
In practice, they carried out explicitly the computation 
of the first four cumulants. 

To draw a connection to their work we note that 
whereas the parametrization x = arctanh-\/i? of the mag- 
nitude of the reficction coefficient is the natural one 
from a geometrical point of view according to Eqs. (|2.9|l . 
and (|?T^ . it is 



1 < li := 2a; - 



\ng 



In 1/ 



(5.8a) 



that behaves statistically as the logarithm of the conduc- 
tance g (or the logarithm of the LDOS i' according to the 
Borland conjecture, see Ref. l53|) in the regime defined by 
the condition 



0<ff < 1. 



(5.8b) 



We denote by U{u, <p', L) the joint probability distribution 
of u and 0, and introduce s = 2L/£. We infer from Eq. 
(|3.15|l that U{u, 0; L) obeys the Fokker-Planck equation 



J 



1 3, ^\dU 

2 2 I ou 



1 d'^U 
-(l-Ccos2</)) — +Csin(20) 



dud(j>^ 



(5.9a) 
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where 



'C cos 20 — 
d 



£ 3, . „ , 
- + -Csin20 



92 
902 



5h 



d 



t:— +Co— + — ( — COS0 ) +— (— sini 



52_ 

5+ 



d 



(5.9b) 



in the large u Umit and in the absence of absorption. 

The fact that can be written as a total derivative plays an essential role when solving, for any given purely 
imaginary z and any given positive integer fc, the eigenvalue equation 



1 



d 



C^-z( -- + -Ccos2</) + Csin(20)— ) + _ (1 - Ccos2(/.) 



fk 



(5.10a) 



for the eigenvalue and the normalized eigenfunction fi^{z, 0) of with periodicity tt that results from insertion 

of the large-deviation ansatz 



U{u,^-L) 



+ CXD 



(5.10b) 



k=0 



into Eq. 

For large u ^ 1, the integration over z is dominated in Eq. 15.10b|l by the region of size ~ 1/m close to the origin 
z = 0. Assume that for a nonvanishing z of order l/u, all eigenvalues iij.{z) are real valued, of descending order, 
and that the largest eigenvalue fJ-j^^giz) is separated from fi/^^iiz) by a finite gap. Consequently, truncation of the 
summation over k to the single term fc = on the right-hand side of Eq. Ij5.10bp produces an exponentially small 
error 



C/(ii,0;L) 



+'°° dz_ r 



(5.11) 



Next, assume the expansions 

+00 +00 

71—1 Tl— 

where the expansion coefficient /Iq"'' encodes the nth cu- 
mulant of u, 

(u")^ = 71! .s + O (s") . (5.13) 

Insertion of the expansion (|5.12() into the Eq. I|5.10a|) for 
fc = can be solved iteratively for n = 0, 1, . . .. For 
ri = 0, Eq. (|5.10a|l reduces to 

= /:^/r, (5.14) 

whose normalized solution is nothing but the stationary 
and normalized solution to Eq. Ij4.4bp . For n = 1, Eq. 
(I5.10a|l reduces to 

A^^'Vr = ^0/o''+(^ - ^Ccos(20) - Csin(2</>)9^) . 

(5.15) 

The expansion coefficient /Xq^'' is obtained from integrat- 
ing Eq. H5.15|l over as a result of being a total deriva- 
tive, 

/x(^) = i(l + c4°0, (5.16a) 



where 

4°^:= r"d(/)COs(20)/^°)(0). (5.16b) 

After comparison of Eqs. H5.13|l and H5.16|l on the one 
hand with Eqs. (|5.3a() and Ij4.15b|l on the other hand, 
we conclude that the first cumulant of the logarithm of 
the conductance is computed exactly to leading order 
in ^/L within the separation-of- variable ansatz for the 
crossover regime (|4.3|) . As a corollary, the crossover (|5.5|) 
for the typical localization length is given exactly by the 
separation-of- variable ansatz. Carrying on the iteration 

to compute the expansion coefficients /Xq"'', however, we 
find that the separation-of-variable ansatz breaks down 
for n = 1,2,... as it does not agree anymore with the 
large-deviation ansatz from Schomerus and Titov in Refs. 
mandlH 



VI. PROBABILITY DISTRIBUTION OF THE 
LDOS 

In this section, the probability distributions P{i'\ L) 
and P{v\ijj) of the LDOS v in the geometries of Figs-Ela. 
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(a) 









ZD 







^ 

y 



L 

(b) 



y 



L 

FIG. 6: (a) Disordered quantum wire of length L closed 
at the right end and connected to a perfect lead at the left 
end. Energy levels within the wire of length L are broadened 
beyond the mean level spacing by the coupling to the perfect 
lead, (b) Disordered quantum wire of length L closed at both 
ends. Energy levels within the wire are broadened beyond the 
mean level spacing by the absorption Im_E = 8" > 0. 



with the boundary conditions associated to the geome- 
tries in Fig. El As long as y is far away from the ends 
of the disordered quantum wire, say in the middle, the 
statistical properties of ^{y) should be independent of 
y. However, the statistical properties of ^{y) depend 
strongly on the boundary conditions at the ends of the 
disordered quantum wire^SiSl For completeness, we will 
understand under the chiral limit of the LDOS H6.1|) 
the property that the Green function (|6.2() anticommutes 
with tT]^, i.e., 

E = go = 9,^0, ^ (6.3) 

whereas we will understand under the standard limit of 
the LDOS liniU the condit ions on the Green function 
(in3 that 

9i = 92 (6-4a) 

or 

.90 > .91 + 92 (6.4b) 



and 03, respectively, are calculated based on the prob- 
ability distribution function of the reflection coefhcient 
r = -v/i? exp(i(/)) obtained in Sec. IIVI 

E:$> 91+92- (6.4c) 



A. Single-particle Green function and reflection 
coefHcient 



For any given realization of the disorder potentials 



12, the LDOS i^iy) is defined bj*^ 
v{y) := --ImtrG(y,y), 



(6.1) 



where the matrix elements G{y, z) in real space of the 
Green function are obtained from solving 



E = £' + i£\ f">0, 
[rj : positive infinitesimal) 



G{y,z) = 5{y - z)cto, 
(6.2) 



Schomerus et al. in Ref. |2l| expressed the probability 
distribution of the LDOS in terms of the probability dis- 
tributions of the reflection coefficient j'l (^r) for outgoing 
plane waves reflected from the region to the left (right) of 
y for the continuous nonrelativistic Schrodinger equation 
in the standard symmetry class (|6.4|) . Here we extend 
their analysis to the continuous relativistic Hamiltonian 
(I2.1a(l . A crucial observation in Ref. I21I is that, for weak 
disorder, the Green function can be expanded in terms of 
free scattering states within a small interval of order of 
the mean free path i since one can assume that there are 
no impurities in this interval. The solution to Eq. 16.2(1 
for y and z belonging to such an impurity-free interval 
can be written as a linear combination of e^^^'+^+^^x+X-j 
gjfc(+y-z)^_^^T^ gZfc(-y+z)^_^T^ e^K-y-z)^_^T^ except 

for the discontinuity at y = z, 



^^-^ky^_^ +^ky ^^+^kz^_^ -^kz 

I 

where x+ '■— (IjO)"^, X- •= (0, l)"^i ©(y) is the Heavi- side function that vanishes when y < and equals one 
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when y > 0, and we have momentarily reinstated the 
Fermi velocity v-p and h in the relativistic dispersion re- 
lation Ke E = £' = hvpk. The relative amplitudes among 
gifc(±i/±2)-^^-^T as the factor (1 — rj^r^)~^ are de- 

termined by taking into account multiple scattering from 
the left and right boundaries of an interval free of impu- 
rities. 

Combining Eqs. H6.2|l and H6.5() reproduces in the rel- 
ativistic limit the mesoscopic relation 

between the LDOS at y and the reflection coefhcients 
and from the disordered region to the left and right 



of y, respectively, derived by Schomerus et al. in Ref. I21I 
for the nonrelativistic continuous Schrodinger equation. 
In the nonrelativistic case, the diagonal matrix element 
of the Green function in real space gives the microscopic 
local density of states. The microscopic local density of 
states exhibits 2fcp oscillations iSi The LDOS (|6.6|l is ob- 
tained after smearing out these oscillations. The effect of 
the relativistic approximation (|2.1a|) is to remove the 2k-p 
oscillations from the outset. Henceforth let us measure 
the LDOS v{y) in units of l/(7rfi,Wp), the DOS of a clean 
wire. 

The probability distribution function P^iy; L^, L^,uj) 
of the LDOS i'{y) defined by Eq. H6.6|l is now simply 
given by 



Jo Jo Jo Jo 



X 5 I 1^ 



1 - Ri^R^ 



1 + i?Li?R - 2^R^R^ cos ( 



, (6.7) 



where (-^r) is the length of the segment of the dis- 
ordered quantum wire to the left (right) of the point y 
where the LDOS is measured. The total length of the 
disordered quantum wire is evidently given by (sec Fig. 

El 



(6.8) 



The dependence of the probability distribution for the 
LDOS on the boundary conditions enters through 
the dependence of the probability distributions 
Pi.{Rh^4>i^'^L^i^)^ and PR(i?R,0R;iR,w) on the 
boundary conditions corresponding to the geometries 
of Fig. QJiSiSi This dependence on boundary conditions 
is implied by the dependence on L^,tjj or i^jW of the 
probability distribution of the two reflection coefficients 
and r^, respectively. When a disordered quantum 
wire of finite length is closed, one must broaden the 
single-particle energy levels beyond the mean level 
spacing to make sense of the energy dependence of the 
mesoscopic LDOS. This is achieved here by introducing 
a finite absorption u = Im E/g_^ > 0, in which case the 
dependence on ^ ^ or Lr 3> i becomes immaterial. 
The absorption uj = ImE/g^ > is switched off when a 
disordered quantum wire of finite length L is connected 
to a perfect lead. We consider first the geometry 
depicted in Fig. 01 and then the geometry depicted in 
Fig.Eb. 



B. Semiopen wire 



To calculate the probability distribution of the LDOS 
in a disordered quantum wire connected to a reservoir 
on the left and closed on the right, we need (i) the joint 
probability distribution ^^(i?^, (f)-^; L^) for the square i?L 
of the magnitude of the reflection coefficient and its 
phase (p-^ on the segment of length to the left of the 
point at which the LDOS is measured and (ii) the proba- 
bility distribution $R(0R;iR) of the phase 0r of the re- 
flection coefficient Tr on the segment of length Lr to the 
right of the point at which the LDOS is measured. The 
calculation of <&r(<^r; ^^r) is outlined in Sec. IIV Al The 
calculation of ^^(i?^, 0^; L^) is outlined in Sec. IIV CI 



1. Chiral and standard symmetry classes 



The probability distribution of the LDOS in the chiral 
symmetry class (|6.3|l is (see the Appendix 



P{v;L^) = Ipyp ( In^zy ) (6.9a) 



8L, 
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The average LDOS is 
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FIG. 7: (a) Three traces for the probability distribution of 
the LDOS are plotted with e = 0, 0.1, +00 (from the bottom 
up at = 0.01) so as to interpolate between the chiral and 
standard symmetry classes 16.31 and 16.41 , respectively. Each 
trace is obtained from numerical integration of Eq. 16.71 for 
the semi-open geometry of Fig.|^ with L-^/l — 2. (b) Three 
traces for the probability distribution of the LDOS are plotted 
with e = 0, 0.1, +00 (from the bottom up at u — 0.1) so as to 
interpolate between the chiral and standard symmetry classes 
16.31 and 16.41 . respectively. Each trace is obtained from 
numerical integration of Eq. 16.71 for the closed geometry of 
Fig. |S|d with Lo = 1/12. We have verified that numerical 
integration of Eq. 16.71 agrees with a representation of Eq. 
16.71 in terms of elementary functions when possible. 



— »2Ll 



(6.10) 



which grows exponentially with Ll- For comparison, 
the probability distribution of the LDOS in the standard 
symmetry class (|6.4f) . which was obtained in Refs.fioland 
I21I is also log-normal. 




(6.11a) 
(6.11b) 

(6.11c) 



Incidentally, the average LDOS is not affected by the 
disorder in the standard symmetry class, (ly) = 1. 
The probability distributions for the LDOS in the chi- 
ral and standard symmetry classes are depicted in Fig. 
HJi. In the thermodynamic limit — > 00, the logarithm 
of the LDOS is self-averaging in the standard symmetry 
class whereas this is not the case in the chiral symmetry 
class. This difference could have been anticipated from 
the identification Inz/ ^ ~2x, valid when L-^/£ 3> 1, with 
the radial coordinate x that obeys the Fokker-Planck 
equations (|2.12|l and H2.13|l . respectively. 



2. Separation ansatz for the crossover regime 

Numerical evaluation of the integrals on the right hand 
side of Eq. (|6.7|l is depicted in Fig. EJi for the crossover 
regime (14 .bf) within our approximation which is encoded 
by the separation-of- variable ansatz ()4.27|l . ()4.28|l . H4.29|l . 
and H4.9|) . 



Closed wire 



To calculate the probability distribution of the LDOS 
This is the probability distribution of a log-normal dis- in a disordered quantum closed at both ends with a finite 
tributed random variable with absorption lu, we need the joint probability distributions 

.00 Pj^{Rj^,(j>j^;L^,uj) and PR (i?R, Lj^, w). We can safely 

(Inz/)^ := / di^P{i^;L^)lnv ignore the dependence on L^. 

in the regime L^^^ 3> £ 

as is done in Sec. HVBl 



0, (6.9b) 



((lni.)2)^ / d,,P{i.;L^) {In ly) 



2 



(6.9c) 



Chiral and standard symmetry classes 
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The probability distribution of the LDOS in the chiral symmetry class (|6.3|) (see appendix 

2 



1 



1 



■ exp 



^(l + i.2)] ^^|1-^2| 

^^0 



2v 



2v 



where K^{x), n — 0,1, . . . are the modified Bessel functions and has the asymptotics 

z^~^/^ exp (— oj/i') , for <C ijj/2, 
P{u;u;) cx I -ln(L..|l-z.|), for |1 - « i 
v^"^/^ exp {—uji') , fori^»2/ijj. 



(6.12a) 



(6.12b) 



For comparison, the probability distribution of the LDOS in the standard symmetry class ()6.4f) . which was obtained 
in Refs.lliandlM 



P(i.;w) 



du- 



exp [—Lu{u — 1)] 



V27ri^3/2 7(i + „.)/(2^) ^u- (l + I/2)/(2i,) 

has the asymptotics 

z^^2 exp (— tj/zy) , for I' ^ ijj/2 



wA'o (^wVu^ - l) + - IK^ (uVu^ - l) , (6.13a) 



P{v;uo) oc 



z^ -'^^2 Qxp (— (jj;y) , for I' 2/(jj. 



(6.13b) 



The probability distributions for the LDOS in the chiral 
and standard symmetry classes are depicted in Fig. Et"- 
Spectral weight from the tails z/ <C ct;/2 and v 2/uj 
of the probability distribution in the standard symmetry 
class is redistributed around z/ = 1 in the chiral symmetry 
class. Consequently, in the limit lo ^ {), the mean value 
of the mesoscopic LDOS in the chiral symmetry class 
becomes 



1 



[Ao(c./2)]2 
"2 



dv 1 



1 



X exp 



1 



(u;/2)[ln(u;/2)]2 



(6.14) 



This result is identical to the smeared Dyson singularity 
of the DOS in the chiral symmetry class, 



Dyso: 



de-. 



2ti 



1 



ln|e/2| 



(c./2)[ln(^/2)]2 



(6.15) 



2. Separation ansatz for the crossover regime 

Numerical evaluation of the integrals on the right hand 
side of Eq. (|6.7|l is depicted in Fig. |7|d for the crossover 
regime ()4.3|l within our approximation which is encoded 
by the separation-of- variable ansatz H4.14|l . (|4.15|) and 
(Oil . 



VII. CONCLUSION 



The probability distribution of the mesoscopic local 
density of states (LDOS) v for a strictly one-dimensional 
problem of Anderson localization with chiral symmetry 
(the chiral symmetry class) was computed in closed form 
for two simply connected geometries assuming a weak 
disorder. As is the case when the chiral symmetry is max- 
imally broken (the standard symmetry class), the prob- 
ability distribution of v strongly depends on boundary 
conditions. In a semiopen geometry we found that the 
probability distribution of v is log-normal with a vanish- 
ing mean as opposed to log-normal with a finite mean 
in the standard symmetry class. In a closed geometry 
with absorption we found that the probability distribu- 
tion of I' has a double-peak structure whereby the sec- 
ond peak turns out to be a logarithmic singularity at 
1/ ^ 1. We verified that the smeared Dyson singularity 
is reproduced by the mean value of the LDOS when the 
absorption is sufficiently small. Furthermore, we found 
one exact and proposed two approximate solutions to 
the functional renormalization group equation obeyed by 
the joint probability distribution for the squared modu- 
lus and the phase of the reflection amplitude of a finite 
and possibly dissipative wire with semiopen and open 
boundary conditions, respectively, from which we could 
extract the approximate crossover of the LDOS v and 
conductance g between the chiral and standard symme- 
try classes. We could show that our approximation is 
exact for the first cumulant but fails to describe higher 
cumulants of the logarithms of v and g in the crossover 
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regime. 
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APPENDIX A: PROBABILITY DISTRIBUTION 
OF THE LDOS IN THE CHIRAL SYMMETRY 
CLASS 

In the appendix we prove Eqs. (|6.9a|) and H6.12a(l . We 
make use of the fact that the probability distributions 
for the squared modulus and phases of the reflection co- 
efhcients r^^^^ = tanhx^^^exp(i(l>j^^^) factorize in the 
chiral symmetry class. We have chosen here to repre- 
sent < (-Rl/r)"'^^^ < 1 by tanhx^/pj. We thus need to 
compute 



poo /"Oc /'2tt /'27r 

■^R' "^Ri -^R' ^ 

Jo Jo Jo Jo 



X ^ — 



1 — tanh^ tanh^ 



1 + tanh tanh — 2 tanhx^ tanhsj^ cos (0^ -|- 0]^) 



where 



^L/r(2;l/R>'/'l/r;-^L/R>^) = '^L/R(a;L/R;^L/R>^) 2 ('?^L/R " o) + (^^L/R " 

Integration over the angles (/>L/p{ yields 



/•OC /"OO 

Jo Jo 



1 

X - 

2 



, , 1 — tanhxr tanh a; r, \ ^ / 1-1- tanhsr tanh a; r, 
6{ ly- ; — ^ — - \+6\v - ^ 



1 -|- tanh tanh x-^ 



1 — tanh x^ tanh x^ 



(Ala) 



(Alb) 



(A2) 



1. Semiopen wire 



2. Closed wire 



Equation (|2.13|) implies that we need to evaluate Eq. 
21 with 



-^Rl^^Ri^R) = (^lajR - oo). 



exp 



(2iL/^) 



,(A3a) 
(A3b) 



It is found that 



1 /"^ 

^(^;iL) = / dxA'L(x;iL)%-exp(2x)). 



(A4) 



Equation (j4.11b|) implies that we need to evaluate Eq. 
with 

A-Jxl;^) = 2ACe-('^/^)™^'^^^S (A5a) 

A'r(xr;c^) = 2ACe-("/')™'^'''"=^, (A5b) 
where = \Kq(u: I2)Y^ . It is found that 



oo /-oo ^— 



(A6) 



Integration over a; yields Eq. H6.9a|l . 



Integration over -p and q yields Eq. (|6.12a|l . 
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